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Abstract 

We calculate the next-to-leading-order cross section for the inclusive production of 
B mesons in pp collisions in the general-mass variable-flavor-number scheme, an 
approach which takes into account the finite mass of the b quarks. We use realistic 
evolved nonperturbative fragmentation functions obtained from fits to e + e~ data 
and compare our results for the transverse-momentum and rapidity distributions at 
a center-of-mass energy of 7 TeV with recent data from the CMS Collaboration at 
the CERN LHC. We find good agreement, in particular at large values of px- 

PACS: 12.38.Bx, 12.39.St, 13.85.Ni, 14.40.Nd 



1 Introduction 



Since the late eighties there has been much interest in the study of 5-meson production 
in pp and pp collisions at hadron colliders, both experimentally and theoretically. The 
first measurements were performed more than two decades ago by the UAl Collaboration 
at the CERN SppS collider pQ operating at a center-of-mass energy of y/S = 0.63 TeV. 
More recent measurements were made by the CDF and DO Collaborations at the Fermilab 
Tevatron running at y^S = 1.8 TeV [21 E] and 1.96 TeV g]. Just recently, the CMS 
Collaboration at the CERN LHC collider published first results for inclusive B + - [5], 
B°- [Bj, and 5 s -meson [7] production in pp collisions at yS = 7 TeV. B + mesons were 
reconstructed via their decay B + — > J/tpK + followed by J/ip — > whereas B° 

mesons were identified through the observation of J / ipK® final states with the subsequent 
decays J/ip — > and K® — > 7r + 7r~. In the case of B s mesons, the reconstructed final 

states were generated by the decay chain B s — > J/ip(fi, J/ij) — > and — > K + K~ . 

From all these measurements the differential cross sections da/dpx and da/dy as well as 
the integrated cross section for p? > 5 GeV (for B + and B° mesons) or p? > 8 GeV (for 
B s mesons) were reported. 

The general-mass variable-flavor-number (GM-VFN) scheme provides a rigorous theo- 
retical framework for the description of the inclusive production of single heavy-flavored 
hadrons, combining the fixed-flavor-number (FFN) [8] and zero-mass variable-flavor-number 
(ZM-VFN) j9] schemes, which are valid in complementary kinematic regions, in a unified 
approach that enjoys the virtues of both schemes and, at the same time, is bare of their 
flaws. Specifically, it resums large logarithms by the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) evolution of nonperturbative fragmentation functions (FFs), guarantees 
the universality of the latter as in the ZM-VFN scheme, and simultaneously retains the 
mass-dependent terms of the FFN scheme without additional theoretical assumptions. It 
was elaborated at next-to-leading order (NLO) for photoproduction [10] and hadropro- 
duction [TT] of charmed hadrons as well as for their production by e + e~ annihilation [T2] . 
It was also applied to obtain predictions for 5-meson hadroproduction [13], which could 
be compared with recent CDF data [1]. An earlier implementation of such an interpo- 
lating scheme is the so-called fixed-order-next-to-leading-logarithm (FONLL) approach, 
in which the conventional cross section in the FFN scheme is linearly combined, with 
the help of a p-r-dependent weight function, with a suitably modified cross section in the 
ZM-VFN scheme implemented with perturbative FFs [T3] . 

In Ref. [13J, nonperturbative FFs for the transitions a — > B, where a is any parton, 
including b and b quarks, were extracted at NLO in the MS factorization scheme with 
rif = 5 flavors from the scaled-energy (x) distributions da/dx of e + e~ — > B + X measured 
by the ALEPH [15] and OPAL [16] Collaborations at the CERN LEP1 collider and by 
the SLD Collaboration [17J at the SLAC SLC collider. As explained in Ref. [13], these 
FFs may be consistently used in our GM-VFN framework. Working at NLO in the GM- 
VFN scheme with these 5-meson FFs, we found excellent agreement with recent CDF 
measurements of da/dpT for pp — > B + X [1], especially in the upper pt range, pr ~ 10 
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GeV [13] . 

The content of this paper is as follows. In Sec. 2, we summarize our input choices of PDFs 
and .B-meson FFs. In Sec. 3, we compare the predictions of the GM-VFN scheme with 
the CMS data from the recent LHC run at \/S = 7 TeV [3 El [7] . Our conclusions are 
given in Sec. 4. 



2 Input PDFs and S-meson FFs 

As PDFs for the proton, we choose one of the most recent parametrizations of the CTEQ 
Collaboration, set CTEQ6.6M [TH], which provides an improvement over the earlier version 
CTEQ6.5M. Both sets were obtained in the framework of a general-mass scheme using 
the input values m c = 1.3 GeV, nib = 4.5 GeV, and a s (mz) = 0.118. In both set, the 
6-quark PDF has its starting scale at /io = rn.&. 

The nonperturbative FFs describing the transition of the b and b quarks into a B me- 
son can be obtained only from experiment. In our earlier work on inclusive 5-meson 
production at the Tevatron [13], we constructed such FFs using as input recent precise 
measurements of the cross section of inclusive S-meson production in e + e~ annihilation 
obtained by the ALEPH |15j . OPAL [IB], and SLD [T7] Collaborations^] These data were 
taken on the Z-boson resonance, so that finite-m;, effects, being of relative order m^/m|, 
are strongly suppressed, which means that we are in the asymptotic regime where the GM- 
VFN scheme is equivalent to the ZM-VFN scheme. The combined fit to the three data sets 
was performed using the NLO value A^=, = 227 MeV corresponding to af^niz) = 0.1181, 
values adopted from Ref. [IB]. The renormalization and factorization scales were chosen 
to be fiR = fiF = Tnz- In accordance with the chosen PDFs, the starting scale of the 
b — >■ B FF was taken to be /io = nib, while the g,q — » B FFs, where q denotes the light 
quarks including the charm quark, were taken to vanish at fiF = /^o- 

For fitting the data, we actually employed two different parametrizations for the b — >■ B 
FF at /io = nib, namely the Peterson ansatz [20J and the simple power ansatz [21]. It 
turned out that the Peterson ansatz led to a very poor fit. Therefore, we shall use in 
this work only the FFs obtained with the power ansatz, whose parameters at the starting 
scale are listed in Table 1 of Ref. [32] . A comparison of the fit performed using this ansatz 
with the three input data sets may be found in Fig. 1 of that reference. 

We note that the data from OPAL and SLD included all 5-hadron final states, in particu- 
lar those with Af, hadrons, while, in the ALEPH analysis, only final states with identified 
B ± and B° mesons were taken into account. Our fit was based on the assumption that 
the FFs of all b hadrons had the same shape. The branching fraction of b — > B + was 
taken equal to that of b — > B° and fixed to 0.397. In our calculations for S s -meson pro- 
duction to be presented below, we shall use the same FFs and rescale them by the factor 

1 Recently, similar data became available also from the DELPHI Collaboration [19] . 
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0.113/0.401, which uses the up-to-date values for the b — > B + and b — >■ B s branching 
fractions quoted by the Particle Data Group [22J. 

We should emphasize that, in the analysis of the available e + e~ annihilation data, the 
charged and neutral B mesons were not separated. Furthermore, the charged states B + 
and B~ could not be distinguished. The FFs obtained in Ref. [T3] are, therefore, valid 
for the average of B + and B~ and, similarly, for the average of B° and B°. 

The factorization scales related to the initial- and final-state singularities entering the 
PDFs and FFs, respectively, can in principle be chosen independently. We checked, 
however, that when estimating theoretical error bands by varying these scales by factors 
of 2 up and down, the extreme values are indeed obtained when the initial- and final-state 
factorization scales are identified. Our default choice of renormalization and factorization 
scales is fifi = fip = rriT = \jp\ + m b- Theoretical uncertainties will be estimated by 
setting fj, R = ^rttit and /ip = ^f^t, and varying £ R and £p about their default values 
£ R = £f = 1 by factors of 2 up and down, restricting the ratio to the range 1/2 < £r/£f < 
2. 

3 Theoretical Predictions for pp — >• B + X and Com- 
parisons with CMS Data 

To obtain an overview of the pt dependence of da/dpx, we first show results for this 
observable, integrated over \y\ < 2.4, for the case of B + production in the GM-VFN 
scheme as described above. This differential cross section is shown in Fig. [I] (left) for 
Pt values between 5 and 30 GeV and in Fig. [I] (right) for larger p T values, up to 100 
GeV, where we expect data to come in the near future when the LHC experiments are 
accumulating more statistics. 

In the pt range between 5 and 30 GeV, the cross section falls off by three orders of 
magnitude. This is essentially due to the behavior of the PDFs as a function of the 
scaling variable x and less so from the behavior of the partonic cross sections. Towards 
low pt values, both the upper edge of the error band and the cross section for the default 
choice of scales rise steadily with decreasing p T value, down to pt = 5 GeV. This is caused 
by the scale dependence of the 6-quark PDF and the FFs. With our choice of scales, they 
fade out and quench the cross section, leading to a turn-over of the pt distributions only 
at pt = and not already at some finite pt value. The lower edge of the error band is 
obtained for £p = 0.5. Here, both the 6-quark PDF and the FFs vanish at pt ~ 8 GeV, 
corresponding to //p = mj = 4.5 GeV. The line representing the lower edge of the error 
band therefore stops at this point. 

The CMS Collaboration measured the differential cross section da/dpT for the production 
of B + mesons [5] (actually the average of B + and B~ mesons), integrated over the y range 
\y\ < 2.4, as a function of pt- The measurement covered the pt range from 5 GeV to 
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30 GeV with five bins. In addition, the differential cross section da/d\y\, integrated over 
the considered range, was given for six \y\ bins. In Ref. [6], the results of the measure- 
ment of _B°-meson production (again for the average of the charge-conjugate states B° 
and B°) were presented. They comprise the differential cross section da/dpT, integrated 
over the y range \y\ < 2.2, in five pr bins between pt = 5 GeV and pr = 40 GeV and 
da/d\y\, integrated over the considered pr range, in five \y\ bins. Since, in this second 
analysis, a larger luminosity was already available, the B° data extend to larger p? values. 

In order to facilitate the comparisons with the CMS measurements [SI E], we integrate 
over the bins using the same binnings. The px bins for B + - and 5°-meson production are 
the same, except for the largest one. Our results are shown in Figs. [2]and[3j where they 
are compared with the experimental data. The errors of the experimental data points are 
obtained from Ref. [6] by adding in quadrature the statistic and systematic errors quoted 
there. The differences between the predictions in Figs. [2] and [3] are entirely due to the 
different bin choices, the FFs being the same in both cases. 




Figure 1: da/dp T [nb/GeV] for pp B + + X at ^/S = 7 TeV in the GM-VFNS. For 
clarity, we split the pt range into a lower part (pt below 30 GeV, left panel) and an 
upper part (pr above 30 GeV, right panel). The central values (solid lines) correspond 
to the default choice of scale parameters, £r = £p = 1. An error band (dashed lines) is 
obtained from variations of the renormalization and factorization scales by factors of 2 up 
and down. The upper end of the error band is reached for £r = 1 and £p = 2 at p? < 21 
GeV and for £r = 0.5 and ^F = latpy>21 GeV, the lower error end is reached for 
£r = 1 and £f = 0.5 at p T < 25 GeV and £ R = 2 and £ F — 1 at p T > 25 GeV. 
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We determine the error band from variations of the scale parameters by factors of 2 as 
described above, except that the minimum of the theoretical prediction is obtained with 
the additional prescription that the FFs are frozen when fip falls below the starting scale 
/i = m b . Otherwise the cross section would become zero for £ F = 0.5 in a large part 
of the first pt bin, so that the lower edge of the error band would become meaningless. 
As is seen in Figs. [2] and |3j the data lie inside the error bands. In the case of B + (B°) 
mesons, the default predictions appreciably overshoot the CMS data in the first three 
(two) pt bins, while they are very close to the CMS data in the residual pt bins. The 
default values of the predicted cross sections are a factor of approximately 2 (1.5) larger 
than the experimental central values in the lowest (next-to- lowest) pt bins. This is caused 
by the fact that, with our choice of scales, large contributions coming from initial-state 
b quarks are present for all finite values of pt- If one changes the factorization scale to 
a lower value, for example by setting £r = 1 and £p = 0.7, the 6-quark PDF vanishes at 
Pt = 4.6 GeV. Furthermore, with our prescription, the PDFs and the FFs are frozen at 
the values they reach at \ip = m;, when pt falls below pt = 4.6 GeV. For this special 
choice of factorization scales, we obtain the cross section values given for the 5°-meson 
case in the column labeled = 1, £p — 0.7 of Tab. 1. For comparison, we present the 
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Figure 2: da/dpp [nb/GeV] (left panel) and da/d\y\ [nb] (right panel) for pp — >■ B + +X at 
NLO in the GM-VFN scheme compared with the CMS data [5]. The central values (solid 
lines) correspond to the choice Cr — — 1- We also show the prediction for da/d\y\ 
obtained with the choice = 1 and £p = 0.7 (dash-dotted line). The error bands (dashed 
lines) are obtained by varying £r and £p by factors of 2 up and down (maximum: £r = 1, 
£f = 2; minimum: £ R — 1, £ F — 0.5). 



6 



da/dpx [nb/GeV] 



pp -> B° + X 
GM-VFNS 
75 = 7 TeV 
M <2.2 
CMS data 



10 15 20 25 



4-10 4 



3-10 4 



2-10 4 



10 4 - 



da/d\y\ [nb] 
GM-VFNS 
p T > 5 GeV 



30 35 40 
Pr [GeV] 



0.5 



pp^ B° + X 

VS = 7 TeV 
CMS data 



1.5 



2 

Ivl 



Figure 3: da/dpx [nb/GeV] (left panel) and da/d\y\ [nb] (right panel) for — >■ £>° + X at 
NLO in the GM-VFN scheme compared with the CMS data [5J. The central values (solid 
lines) correspond to the choice £r = £f = 1- We also show the prediction for da/d\y\ 
obtained with the choice ^r = 1 and £p = 0.7 (dash-dotted line). The error bands (dashed 
lines) are obtained by varying £r and £p by factors of 2 up and down (maximum: £r = 1, 
^ = 2; minimum: = 1, £p = 0.5). 



experimental results in the second column of this table and the default-scale results of 
Fig. [3] (left) in the third one. We see that the theoretical values of the cross sections in the 
five pt bins agree with the experimental values quite well, within the errors. The total 
5°-meson production cross section determined by CMS in the considered kinematic range 
is cr t ot = 33.2 ±4.3 fib. For the default choice of scales = £f = 1, we find a to t = 61.7 jub, 
while the result for £r = 1 and £p = 0.7 is 35.0 fib, in very good agreement with the 
data. A similar comparison may be performed for pp — > B + X, with similar conclusions, 
as can be inferred from Fig. [2] (right panel), where we show the corresponding results for 
da/d\y\. The theoretical predictions are almost identical, since the FFs for b — > B + and 
b — >■ B° are taken to be the same and there is only a tiny difference due to the different 
upper ends of the pt ranges. 

As explained above, massless contributions, in particular the ones due to incoming b 
quarks, dominate the total cross section towards low px values. These contributions 
lead to an increase of da/dpx in the limit px — > because the heavy-quark PDFs carry 
resummed logarithms, which are not fully cancelled by the subtraction terms in the GM- 
VFN approach, which are implemented at NLO, i.e. at fixed order only. This increase can 
be tamed by imposing the kinematic cut s > 4ml f° r the partonic center-of-mass energy 
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p T (in GeV) Data [6J £ R = £ F = 1 £r = 1, = 0.7 £ a = 0.2 
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535 ± 68 
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17-24 


145 ± 20 
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141 
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24-40 


27 ±4 


25 


22 
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Table 1: Predictions for the differential cross section da/dpr [nb/GeV] of 5°-meson 
production with different renormalization and factorization scales compared with the CMS 
data [6], for which the statistical and systematic errors are added in quadrature. The 
values presented in the second and third columns are also displayed in Fig. [3] (left). 



s also for the massless contributions. Furthermore, a judicious choice of the factorization 
scale, e.g. 



liF = y™>i + &j%, (1) 

with a parameter £ a < 1, can boost the transition fip — > //o = for pr — > 0. This 
prescription creates a turn-over of the pt distribution towards low pt values and also 
allows us to obtain a reasonable description of the CDF data [3], which were taken at 
lower pt values. The CMS data start at px = 5 GeV, and a turn-over is not visible in 
da/dpT- However, the ansatz of Eq. leads to a reduction of the px distribution for small 
Pt values, i.e. to a significant change of da/dpx in the first two pt bins. The cross section 
values obtained for B° mesons using the scale choice of Eq. ({!]) with £ a = 0.2 are presented 
in the last column of Tab. 1. We find that this approach leads to a better description of 
the CMS data, which is, however, not as good as for the scale choice £f = 0.7 (fourth 
column of Tab. 1). 

As a side remark, we note that the behavior towards small pt values is not due to a 
shift in the average 5-meson to 6-quark momentum fraction. This may be observed by 
calculating the quantity 

Jdzzda(p T ) . , 

{Z)M = fdzda(p T ) > (2) 

where z is the scaling variable of the FFs and it is understood that the integration is also 
done over the rapidity interval \y\ < 2.4 relevant for the CMS measurement [5]. We find a 
rather weak dependence on pr- In fact, (z) decreases from 0.770 at pt = 5 GeV to 0.749 
at pt = 30 GeV, which means that, in our applications, the b — )■ B FF is always probed 
around its maximum (see Ref. [T3]). 

We now discuss the \y\ distributions da/d\y\ of B + and B° production shown in the right 
panels of Figs. [2] and [3j respectively. The bulk of these cross sections comes from the 
lowest pt bin, where the theoretical uncertainties are largest, as is evident from Tab. 
1. However, it is interesting to find out how much the shapes of these differential cross 
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Figure 4: Bda/dpT [nb/GeV] (left panel) and Bda/d\y\ [nb] (right panel) for — > -B s +X 
at NLO in the GM-VFN scheme compared with the CMS data [7J. The branching fraction 
of the decay B s — > J/tfxfr is assumed to be B = 1.3 x 10 -3 The central values (solid 
lines) correspond to the choice £r = £f = 1- The error bands (dashed lines) are obtained 
by varying and £p by factors of 2 up and down (maximum: = 0.5, £p = 1; minimum: 
& = = 0.5). 



sections depend on the various scale choices. In order to get some idea about this, we 
include in the right panels of Figs. [2] and [3] as dot-dashed histograms also the predictions 
evaluated using the scale choice £r = 1 and £p = 0.7, as in the fourth column in Tab. 1. 
They agree fairly well with the CMS data, while the default predictions (£r = £p = 1), 
shown as solid histograms, significantly overshoot the CMS data as expected, but their 
shapes are still reasonable. 

Finally, in Fig. |4j we present our predictions for the production of B s mesons and compare 
them with the experimental data published by the CMS Collaboration in Ref. [7j. da/dpr 
was measured in four pt bins between px = 8 and 50 GeV and integrated over \y\ < 2.4, 
and da/d\y\ was measured in four \y\ bins spanning this \y\ range and integrated over the 
full pt range considered. Both the experimental data and our theoretical predictions refer 
to the product of cross section times branching fraction B for B s — > J / ip<j), for which we 
adopt the value 1.3 x 10~ 3 from Ref. [22]. In this case, we find better agreement between 
theory and experiment over the full pt range, probably due to the fact that very low 
values of p?, with p T < 8 GeV, are excluded from this analysis. The total cross section 
times branching fraction measured by CMS for 8 GeV < px < 50 GeV and \y\ < 2.4 is 
6.9 ± 0.8 nb, while our calculation yields 7.2 nb. 
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4 Conclusions 



In summary, we applied the GM-VFN scheme to obtain NLO predictions for the produc- 
tion of B mesons in pp collisions at the LHC The comparison with experimental data 
from the CMS Collaboration at v^S* = 7 TeV generally shows good agreement between 
theory and experiment, in particular at large pt values. The agreement is particularly 
good for the case of 5 s -meson production, where data are restricted to pt values above 8 
GeV. At low pt values, we observe large scale uncertainties. 

Future data collection at the LHC will allow us to extend the comparisons with theoretical 
predictions to much wider pt ranges. If also the systematic uncertainties can be further 
reduced, we may expect that 5-meson production will play an increasingly important role 
in constraining size and shape of both PDFs and FFs. 
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